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Abstract. We unveil an alluring alternative to parametric search that applies to both 
the non-geodesic and geodesic Frechet optimization problems. This randomized approach 
is based on a variant of red-blue intersections and is appealing due to its elegance and 
practical efficiency when compared to parametric search. 

We present the first algorithm for the geodesic Frechet distance between two polygonal 
curves A and B inside a simple bounding polygon P. The geodesic Frechet decision 
problem is solved almost as fast as its non-geodesic sibling and requires 0{N 2 log k) time 
and 0(k + N) space after O(k) preprocessing, where N is the larger of the complexities 
of A and B and k is the complexity of P. The geodesic Frechet optimization problem is 
solved by a randomized approach in 0(k + N 2 log kN log N) expected time and 0(k + N 2 ) 
space. This runtime is only a logarithmic factor larger than the standard non-geodesic 
Frechet algorithm 0J. Results are also presented for the geodesic Frechet distance in a 
polygonal domain with obstacles and the geodesic Hausdorff distance for sets of points or 
sets of line segments inside a simple polygon P. 



1. Introduction 

The comparison of geometric shapes is essential in various applications including com- 
puter vision, computer aided design, robotics, medical imaging, and drug design. The 
Frechet distance is a similarity metric for continuous shapes such as curves or surfaces 
which is defined using reparametrizations of the shapes. Since it takes the continuity of the 
shapes into account, it is generally a more appropriate distance measure than the often used 
Hausdorff distance. The Frechet distance for curves is commonly illustrated by a person 
walking a dog on a leash pjj. The person walks forward on one curve, and the dog walks 
forward on the other curve. As the person and dog move along their respective curves, a 
leash is maintained to keep track of the separation between them. The Frechet distance is 
the length of the shortest leash that makes it possible for the person and dog to walk from 
beginning to end on their respective curves without breaking the leash. See section [2] for a 
formal definition of the Frechet distance. 
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Most previous work assumes an obstacle-free environment where the leash connecting 
the person to the dog has its length defined by an L p metric. In [4] the Frechet distance 
between polygonal curves A and B is computed in arbitrary dimensions for obstacle-free 
environments in 0{N 2 \ogN) time, where N is the larger of the complexities of A and B. 
Rote [23J computes the Frechet distance between piecewise smooth curves. Buchin et al. [7] 
show how to compute the Frechet distance between two simple polygons. Frechet distance 
has also been used successfully in the practical realm of map matching [26]. All these works 
assume a leash length that is defined by an L p metric. 

This paper's contribution is to measure the leash length by its geodesic distance inside 
a simple polygon P (instead of by its L p distance). To our knowledge, there are only two 
other works that employ such a leash. One is a workshop article [18] that computes the 
Frechet distance for polygonal curves A and B on the surface of a convex polyhedron in 
0(N 3 k 4 log(fciV)) time. The other paper [12] applies the Frechet distance to morphing by 
considering the polygonal curves A and B to be obstacles that the leash must go around. 
Their method works in 0(N 2 log 2 N) time but only applies when A and B both lie on the 
boundary of a simple polygon. Our work can handle both this case and more general cases. 
We consider a simple polygon P to be the only obstacle and the curves, which may intersect 
each other or self-intersect, both lie inside P. 

A core insight of this paper is that the free space in a geodesic cell (see section [2j) is 
x-monotone, y-monotone, and connected. We show how to quickly compute a cell boundary 
and how to propagate reachability through a cell in constant time. This is sufficient to 
solve the geodesic Frechet decision problem. To solve the geodesic Frechet optimization 
problem, we replace the standard parametric search approach by a novel and asymptotically 
faster (in the expected case) randomized algorithm that is based on red-blue intersection 
counting. We show that the geodesic Frechet distance between two polygonal curves inside 
a simple bounding polygon can be computed in 0(k + iV 2 log kN log N) expected time and 
0(k + iV 3 log kN) worst-case time, where N is the larger of the complexities of A and B 
and k is the complexity of the simple polygon. The expected runtime is almost a quadratic 
factor in k faster than the straightforward approach, similar to [12] . of partitioning each 
cell into 0(k 2 ) subcells. Briefly, these subcells are simple combinatorial regions based on 
pairs of hourglass intervals. It is notable that the randomized algorithm also applies to 
the non-geodesic Frechet distance in arbitrary dimensions. We also present algorithms to 
compute the geodesic Frechet distance in a polygonal domain with obstacles and the geodesic 
Hausdorff distance for sets of points or sets of line segments inside a simple polygon. 

2. Preliminaries 

Let k be the complexity of a simple polygon P that contains polygonal curves A and 
B in its interior. In general, a geodesic is a path that avoids all obstacles and cannot be 
shortened by slight perturbations [20] . However, a geodesic inside a simple polygon is simply 
a unique shortest path between two points. Let n(a, b) denote the geodesic inside P between 
points a and b. The geodesic distance d(a, b) is the length of a shortest path between a and 
b that avoids all obstacles, where length is measured by L2 distance. 

Let I, t, and if denote decreasing, increasing, and decreasing then increasing functions, 
respectively. For example, "H is ||-bitonic" means that H is a function that decreases 
monotonically then increases monotonically. A bitonic function has at most one change in 
monotonicity. 
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Figure 1: Shortest paths in the hourglass define H-^ ^. 

The Frechet distance for two curves A, B : [0, f] — > M 1 is defined as 

S F (A,B) = inf sup d'( A(f(t)),B(g(t)) ) 

/,p:[0,l]-»[0,l] te[0,l] 

where / and g range over continuous non-decreasing reparametrizations and d! is a distance 
metric for points, usually the Li distance, and in our setting the geodesic distance. For a 
given e > the free space is defined as FS £ (A,B) = {(s,t) \ d'(A(s), B(t)) < e}£[0, l] 2 . A 
free space cell C C [0, l] 2 is the parameter space defined by two line segments ab G A and 
cd € -B, and the free space inside the cell is FS E (ab,cd) = FS e (A,B) n C. 

The decision problem to check whether the Frechet distance is at most a given e > is 
solved by Alt and Godau [4] using a free space diagram which consists of all free space cells 
for all pairs of line segments of A and B. Their dynamic programming algorithm checks 
for the existence of a monotone path in the free space from (0,0) to (1, f) by propagating 
reachability information cell by cell through the free space. 

2.1. Funnels and Hourglasses 

Geodesies in a free space cell C can be described by either the funnel or hourglass 
structure of |Hj. A funnel describes all shortest paths between a point and a line segment, 
so it represents a horizontal (or vertical) line segment in C. An hourglass describes all 
shortest paths between two line segments and represents all distances in C. 

The funnel T ^ describes all shortest paths between an apex point p and a line segment 

cd. The boundary of J 7 ^ is the union of the line segment cd and the shortest path chains 
ir(p, c) and n(p,d). The hourglass W^g^ describes all shortest paths between two line 
segments ab and cd. The boundary of 7"^;^ is composed of the two line segments ab, 

cd and at most four shortest path chains involving a, 6, c, and d. See Figured! Funnel and 
hourglass boundaries have O(k) complexity because shortest paths inside a simple polygon 
P are acyclic, polygonal, and only have corners at vertices of P |15j . 

Any horizontal or vertical line segment in a geodesic free space cell is associated with 
a funnel's distance function F p ^ : [c, d] — > R with F p ~^{q) = d(p,q). The below three 
results are generalizations of Euclidean properties and are omitted. See [10] for details. 

Lemma 2.1. F p ^ is [^-bitonic. 

Corollary 2.2. Any horizontal (or vertical) line segment in a free space cell has at most 
one connected set of free space values. 

Consider the hourglass TO^ ^ in Figure [H Let the shortest distance from a to any point 

on cd occur at M a G cd. Define similarly. As p varies from a to 6, the minimum distance 
from p to cd traces out a function H-^ ^ : [a, b] — » R with H-^ -^{p) = min gg j C (i ] d(p, q). 
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Lemma 2.3. H-^ ^ is []-bitonic. 

3. Geodesic Cell Properties 

Consider a geodesic free space cell C for polygonal curves A and B inside a simple 
polygon. Let ab G A and cd G B be the two line segments denning C. 

Lemma 3.1. For any e, cell C contains at most one free space region R, and R is x- 
monotone, y -monotone, and connected. 

Proof. The monotonicity of R follows from Corollary 12,21 For connectedness, choose any two 
free space points (pi, q±), (p2, (72), and construct a path connecting them in the free space 
as follows: move vertically from (pi,qi) to the minimum point on its vertical. Do the same 
for (p2,Q2)- By Lemma |2~T1 this movement causes the distance to decrease monotonically. 
By Lemma l2~3l any two minimum points are connected by a jf-bitonic distance function 
H-aS cd ( c ^- section |2~T|) . but as the starting points are in the free space - and therefore have 
distance at most e - all points on this constructed path lie in the free space. ■ 

Given C's boundaries, it is possible to propagate reachability information (see section 
[2]) through C in constant time. This follows from the monotonicity and connectedness of 
the free space in C and is useful for solving the geodesic decision problem. 

4. Red-Blue Intersections 

This section shows how to efficiently count and report a certain type of red-blue inter- 
sections in the plane. This problem is interesting both from theoretical and applied stances 
and will prove useful in section f5T3l for the Frechet optimization problem. 

Let R be a set of m "red" curves in the plane such that every red curve is continuous, 
x-monotone, and monotone decreasing. Let B be a set of n "blue" curves in the plane where 
each blue curve is continuous, x-monotone, and monotone increasing. Assume that the 
curves are defined in the slab [a, (3] x R, and let I{k) be the time to find the at most one 
intersection of any red and blue curveQ 

Theorem 4.1. The number of red-blue intersections between R and B in the slab [a, f3] x R 
can be counted in 0{N log N) total time, where N = max(m,n). These intersections can be 
reported in 0(N log N + K ■ 1(h)) total time, where K is the total number of intersections 
reported. After 0{N log N) preprocessing time, a random red-blue intersection in [a, [3] x R 
can be returned in 0(logN + I{k)) time, and the red curve involved in the most red-blue 
intersections can be returned in 0(1) time. All operations require 0(N) spaced 

Proof Sketch. Figure [2] illustrates the key idea. Suppose a red curve r^(x) lies above a blue 
curve b2(x) at x = a. If it is also true that r^{x) lies below b2(x) at x = (3, then these 
monotone curves must intersect in [a, [3] x R. Two sorted lists L a , Lp of curve values store 
how many blue curves lie below each red curve at x = a and x = (3. Subtracting the values 
in L a and Lp yields the number of actual intersections for each red curve in [a, (3] x R (and 

1 There is at most one intersection due to the monotonicities of the red and blue curves. 

2 Palazzi and Snoeyink [21] also count and report red-blue intersections using a slab-based approach. 
However, their work is for line segments instead of curves, and they require that all red segments are disjoint 
and all blue segments are disjoint. We have no such disjointness requirement. 



GEODESIC FRECHET DISTANCE INSIDE A SIMPLE POLYGON 



197 




Figure 2: r^(x) lies above two blue curves at x = a but only lies above one blue curve at 
x = (3. Subtraction reveals that r^(x) has one intersection in the slab [a,0\ x E. 



also reveals the red curve that is involved in the most intersections). Intersection counting 
simply sums up these values. Intersection reporting builds a balanced tree from L a and Lp. 

To find a random red-blue intersection in [a, ff\ x M, precompute the number k of red- 
blue intersections in [a, /3] x R, Pick a random integer between 1 and k and use the number 
of intersections stored for each red curve to locate the particular red curve r»(ac) that is 
involved in the randomly selected intersection. By searching a persistent version of the 
reporting structure [24], rj(x)'s jth red-blue intersection can be returned in 0(log N + I(k)) 
query time after 0(N log N) preprocessing time. ■ 



5. Geodesic Frechet Algorithm 

5.1. Computing One Cell's Boundaries in O(logfc) Time 

A boundary of a free space cell is a horizontal (or vertical) line segment. This boundary 
can be associated with a funnel T that has a jf-bitonic distance function F p ^ (cf. 
Lemma [2.ip . Given e > 0, computing the free space on a cell boundary requires finding the 
(at most two) values t\ t t% such that F p ^(ti) = F ^(£2) = e (see Figure [3l). 

d 

h spac i 



c 

a) b) c) 

Figure 3: a & b) A funnel T ^ is associated with a cell boundary and has a bitonic dis- 
tance function F p ^. c) The (at most two) values t±, t2 such that F p ^(ti) = 
F ^(£2) = e define the free space on a cell boundary. 




Lemma 5.1. Both the minimum value of F p ^ and the (at most two) values t\, £2 such that 
F p Bd(^l) = crf(* 2 ) = £ can ^ e f oun d for any e > in 0(\ogk) time (after preprocessing). 

Proof Sketch. After 0{k) shortest path preprocessing [HUH], a binary search is performed 
on the 0(k) arcs of F p ^ in O(logZc) time. See our full paper [10] for details. ■ 
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Figure 4: Critical values of the Frechet distance 



Corollary 5.2. The free space on all four boundaries of a free space cell can be found in 
0(log k) time by computing t\ and t2 for each boundary. 

5.2. Geodesic Frechet Decision Problem 

Theorem 5.3. After preprocessing a simple polygon P for shortest path queries in 0{k) 
time [13J ; the geodesic Frechet decision problem for polygonal curves A and B inside P can 
be solved for any e > in 0(N 2 log A;) time and 0(k + N) space. 

Proof. Following the standard dynamic programming approach of [4], compute all cell bound- 
aries in 0{N 2 log k) time (cf. Corollary 15. 2\ . and propagate reachability information through 
all cells in 0(N 2 ) time. 0(k) space is needed for the preprocessing structures of [13J, and 
only 0(N) space is needed for dynamic programming if two rows of the free space diagram 
are stored at a time. ■ 

5.3. Geodesic Frechet Optimization Problem 

Let e* be the minimum value of e such that the Frechet decision problem returns true. 
That is, e* equals the Frechet distance 5f(A, B). Parametric search is a technique commonly 
used to find e* (see PI BUSIES] )H The typical approach to find e* is to sort all the cell 
boundary functions based on the unknown parameter e*. The comparisons performed during 
the sort guarantee that the result of the decision problem is known for all "critical values" 
[4] that could potentially define e*. Traditionally, such a sort operates on cell boundaries 
of constant complexity. The geodesic case is different because each cell boundary has O(k) 
complexity. As a result, a straightforward parametric search based on sorting these values 
would require 0(kN 2 log kN) time even when using Cole's [9] optimization^ 

We present a randomized algorithm with expected runtime 0(k + N 2 log kN log N) and 
worst-case runtime 0(k + N 3 log kN). This algorithm is an order of magnitude faster than 
parametric search in the expected case. 

Each cell boundary has at most one free space interval (cf. Lemma [2TT]) . The upper 
boundary of this interval is a function bij(e), and the lower boundary of this interval is a 
function a^ie). See Figure UK The seminal work of Alt and Godau [4] defines three types 

^An easier to implement alternative to parametric search is to run the decision problem once for every 
bit of accuracy that is desired. This approach runs in 0(BN 2 log k) time and 0(k + N) space, where B is 
the desired number of bits of accuracy |25| . 

4 A variation of the general sorting problem called the "nuts and bolts" problem (see |17| 1 is tantalizingly 
close to an acceptable 0(N 2 log N) sort but does not apply to our setting. 
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of critical values that are useful for computing the exact geodesic Frechet distance. There 
are exactly two type (a) critical values associated with distances between the starting points 
of A and B and the ending points of A and B. Type (b) critical values occur 0{N 2 ) times 
when aij(e) = bij{e). See Figure [4b. Type (a) and (b) critical values occur 0(N 2 ) times 
and are easily handled in 0(N 2 log k log N) time. This process involves computing values in 
0(N 2 log/c) time, sorting in 0(N 2 log N) time, and running the decision problem in binary 
search fashion 0(log N) times. Resolving the type (a) and (b) critical values as a first step 
will simplify the randomized algorithm for the type (c) critical values. 

Alt and Godau [4] show that type (c) critical values occur when the position of aij{e) 
in cell Cij equals the position of bkj{e) in cell Ckj in the free space diagram. See Figure 
Hh. As e increases, by Lemma l2~Tl ctij(s) is j-monotone on the cell boundary and 6y(e) is 
J-monotone (see Figure |4)d). As illustrated in Figured];, o>ij(e) and bkj(e) intersect at most 
once. This follows from the monotonicities of aij(e) and bkj(e). Hence, there are 0(N 2 ) 
intersections of a,ij(s) and &jy(e) in row j and a total of 0(N 3 ) type (c) critical values over 
all rows. There are also 0(N 2 ) intersections of ay(e) and 6jfc(e) in column i and a total of 
0(iV 3 ) additional type (c) critical values over all columns. 

Lemma 5.4. The intersection of oy(e) and bfcj(e) can be found for any e > in 0(logk) 
time after preprocessing. 

Proof Sketch. Build binary search trees for a,ij(s) and bki(e) and perform a binary search. 



Theorem 14.11 requires that all ajj(e) and b^i{e) are defined in the slab [a, (3] x M that 
contains e*. Precomputing the type (a) and type (b) critical values of [4] shrinks the slab 
such that no left endpoint of any relevant a>ij(e), bki(e) appears in [a, 0\ X K when processing 
the type (c) critical values. In addition, (Hj{e), bki{e) can be extended horizontally so that no 
right endpoint appears in [a, (3] x M. These changes do not affect the asymptotic number of 
intersections and allow Theorem 14. II to count and report type (c) critical values in [a,/3] xR. 

The below randomized algorithm solves the geodesic Frechet optimization problem in 
0(k + N 2 log kN log N) expected time. This is faster than the standard parametric search 
approach which requires 0(kN 2 log kN) time. 

Randomized Optimization Algorithm 

(1) Precompute and sort all type (a) and type (b) critical values in 0{N 2 \ogkN) time 
(cf. Lemma l57i]) . Run the decision problem OilogN) times to resolve these values 
and shrink the potential slab for e* down to [a,/3] xRin 0(N 2 log k log N) time. 

(2) Count the number Kj of type (c) critical values for each row j in the slab [a, 0\ x M. 
using Theorem 14. 1[ Let Cj be the resulting counting data structure for row j. 

(3) To achieve a fast expected runtime, pick a random intersection "dj for each row using 
CjE See Theorem ED 

(4) To achieve a fast worst-case runtime, use Cj to find the a^-(e) curve in each row that 
has the most intersections (see Theorem 14.11) . Add all intersections in [a, /9] xl that 
involve qmj(s) to a global pool V of unresolved critical valued and delete aAfj(e) 
from any future consideration. 

^Picking a critical value at random is related to the distance selection problem [6] and is mentioned in [2], 
but to our knowledge, this alternative to parametric search has never been applied to the Prechet distance. 
6 The idea of a global pool is similar to Cole's optimization for parametric search [9]. 
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(5) Find the median H of the values in V in 0(N 2 ) time using the standard median 
algorithm mentioned in [17]. Also find the median of the O(N) randomly selected 
"&j in O(N) time using a weighted median algorithm based on the number of critical 
values kj for each row j. 

(6) Run the decision problem twice: once on H and once on ^. This shrinks the search 
slab [a,0] x R and at least halves the size of V. Repeat steps [2] through [6] until all 
row-based type (c) critical values have been resolved. 

(7) Resolve all co/i«mn-based type (c) critical values in the same spirit as steps [2] through 
[6] and return the smallest critical value that satisfied the decision problem as the value 
of the geodesic Frechet distance. 

Theorem 5.5. The exact geodesic Frechet distance between two polygonal curves A and B 
inside a simple bounding polygon P can be computed in 0{k + N 2 log kN log N) expected 
time and 0(k + N 3 logkN) worst-case time, where N is the larger of the complexities of A 
and B and k is the complexity of P. 0(k + N 2 ) space is required. 

Proof. Preprocess P once for shortest path queries in 0(k) time [13]. In the expected case, 
each execution of the decision problem will eliminate a constant fraction of the remaining 
type (c) critical values due to the proof of Quicksort's expected runtime and the median of 
medians approach for ty. Consequently, the expected number of iterations of the algorithm 
is 0(logiV 3 ) = O(logiV). 

In the worst-case, each of the O(N) aij(e) in a row will be picked as aMj(e)- Therefore, 
each row can require at most O(N) iterations. Since all rows are processed each iteration, 
the entire algorithm requires at most O(N) iterations for row-based critical values. By a 
similar argument, co/wmn-based critical values also require at most O(N) iterations. 

The size of the pool V is expressed by the inequality S(x) < 5 '( a ~ 1 )+ (:) ( JV ) ; w here x 
is the current step number, and 5(0) = 0. Intuitively, each step adds 0(N 2 ) values to V 
and then at least half of the values in V are always resolved using the median S. It is not 
difficult to show that S(x) G 0(N 2 ) for any step number x. 

Each iteration of the algorithm requires intersection counting and intersection calcula- 
tions for O(N) rows (or columns) at a cost of 0(N 2 log kN) time. In addition, the global 
pool V has its median calculated in 0(N 2 ) time, and the decision problem is executed in 
0(N 2 log A;) time. Consequently, the expected runtime is 0(k + N 2 log kN log N) and the 
worst-case runtime is 0(k + N 3 log kN) including 0(k) preprocessing time [13] for geodesies. 
The preprocessing structures use 0{k) space that must remain allocated throughout the al- 
gorithm, and the pool V uses 0(N 2 ) additional space. ■ 

Although the exact non-geodesic Frechet distance is normally found in 0(N 2 log N) time 
using parametric search (see [4]), parametric search is often regarded as impractical because 
it is difficult to implement^ and involves enormous constant factors [9]. To the best of our 
knowledge, the randomized algorithm in section 15.31 provides the first practical alternative 
to parametric search for solving the exact non-geodesic Frechet optimization problem in R'. 

Theorem 5.6. The exact non-geodesic Frechet distance between two polygonal curves A 
and B in M. 1 can be computed in 0(N 2 log 2 N) expected time, where N is the larger of the 
complexities of A and B. 0(N 2 ) space is required. 



Quicksort-based parametric search has been implemented by van Oostrum and Veltkamp [25j using a 
complex framework. 
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a) b) c) 



Figure 5: a) A funnel for a <5c-cell can be found by extending a cell's initial leash along one 
segment to create a path sketch and then b) snapping this sketch into a homotopic 
shortest path, c) A funnel T Q ^ has O(kN) complexity, but the distance function 
F Q ^ has only 0(k) complexity because d(o,p) is a constant. 

Proof. The argument is very similar to the proof of Theorem 15.51 The main difference is 
that non-geodesic distances can be computed in 0(1) time (instead of O(logfc) time). ■ 

6. Geodesic Frechet Distance in a Polygonal Domain with Obstacles 

Consider the real-life situation of a person walking a dog in a park. If the person and 
dog walk on opposite sides of a group of trees, then the leash must go around the trees. More 
formally, suppose the two polygonal curves A and B lie in a planar polygonal domain V [19] 
of complexity k. The leash is required to change continuously, i.e., it must stay inside T> and 
may not pass through or jump over an obstacle. It may, however, cross itself. Let 5c be the 
geodesic Frechet distance for this scenario when the leash length is measured geodesicallyll 

Due to the continuity of the leash's motion, the free space inside a geodesic cell is 
represented by an hourglass - just as it was for the geodesic Frechet distance inside a simple 
polygon. Hence, free space in a cell is x-monotone, y-monotone, and connected (cf. Lemma 
13.1(1 . and reachability information can be propagated through a cell in constant time. 

The main task in computing 5c is to construct all cell boundaries. Once the cell bound- 
aries are known, the decision and optimization problems can be solved by the algorithms 
for the geodesic Frechet distance inside a simple polygon (cf. Theorems 15.31 and 15. 5|) . We 
use Hershberger and Snoeyink's homotopic shortest paths algorithm [16j to incrementally 
construct all cell boundary funnels needed to compute 5c- To use the homotopic algorithm, 
the polygonal domain V should be triangulated in O(klogk) time [19], and all obstacles 
should be replaced by their vertices. A shortest path map [19] can find an initial geodesic 
leash Lj between the start points of the polygonal curves A and B in 0(k\ogk) time. 

Lemma 6.1. Given the initial leash for the bottom-left corner of a 5c~cell C, all four funnel 
boundaries of C and the initial leashes for cells adjacent to C can be computed in O(k) time. 

Proof. The funnels representing cell boundaries are constructed incrementally. The idea is 
to extend the initial leash into a homotopic "sketch" that describes how the shortest path 
should wind through the obstacles and then to "snap" this sketch into a shortest path (see 
Figures [5h and [5b). 



'We recently learned that this topic has been independently explored in [8]. 
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Homotopic shortest paths have increased complexity over normal shortest paths because 
they can loop around obstacles. For example, if the person walks in a triangular path 
around all the obstacles, then the leash follows a homotopic shortest path that can have 
O(k) complexity in a single cycle around the obstacles. By repeatedly winding around the 
obstacles O(N) times, a path achieves O(kN) complexity. 

To avoid spending O(kN) time per cell, we extend a previous homotopic shortest path 
into a sketch by appending a single line segment to the previous path (see Figured). Adding 
this single segment can unwind at most one loop over a subset of obstacles, so only the most 
recent 0(k) vertices of the sketch will need to be updated when the sketch is snapped into 
the true homotopic shortest path. A turning angle is used to identify these O(k) vertices by 
backtracking on the sketch until the angle is at least 2ir different from the final angle. 

Putting all this together, a boundary for a free space cell can be computed in O(k) 
time by starting with an initial leash Lj of O(kN) complexity, constructing a homotopic 
sketch by appending a single segment to Lf, backtracking with a turning angle to find O(k) 
vertices that are eligible to be changed, and finally "snapping" these 0(k) vertices to the 
true homotopic shortest path using Hershberger and Snoeyink's algorithm [16]. The result 
is a funnel that describes one cell boundary. 

By extending Lj in four combinatorially distinct ways, all four cell boundaries can be 
defined. Specifically, we can extend Lj along the current ab G A segment to form the first 
funnel or along the cd G B segment to form the second funnel. The third funnel is created 
by extending Lj along ab G A and then cd G B. The fourth funnel is created by extending 
Lj along cd G B and then ab G A. These cell boundaries conveniently define the initial leash 
for cells that are adjacent to C. m 

Theorem 6.2. The 5c decision problem can be solved in 0(kN 2 ) time and 0(k + N) space. 

Proof. Each cell boundary is a funnel J 7 ^ with 0(kN) complexity [H]. However, this high 
complexity is a result of looping over obstacles, and most of these points do not affect the 
funnel's distance function F Q -g. As illustrated in Figure [5b, F o ^ has only 0(k) complexity 
because only vertices 7r(p, c) U ir(p, d) contribute arcs to F o ^. 

Construct all cell boundary funnels in 0(kN 2 ) time (cf. Lemma l6T|) . intersect each 
funnel's distance function with y = e in 0(N 2 logk) time, and propagate reachability in- 
formation in 0{N 2 ) time. Only 0(k + N) space is needed for dynamic programming when 
storing only two rows at a time. ■ 

Theorem 6.3. The 5c optimization problem can be solved in 0(kN 2 + N 2 log kN log N) 
expected time and 0(kN 2 ) spaced 

Proof. The 5c optimization problem can be solved using red-blue intersections. C^logA^) 
steps are performed in the expected case by Theorem 15.51 Each step has to perform in- 
tersection counting in 0(N 2 log kN) time and solve the decision problem. If the funnels 
are precomputed in 0(kN 2 ) time and space, then the decision problem can be solved in 
0(N 2 logk) time. Hence, after 0(kN 2 ) time and space preprocessing, 5c can be found in 
0(logN) expected steps where each step takes 0(N 2 log kN) time. ■ 

9 If space is at a premium, the algorithm can also run with 0(k + N 2 ) space and 0(kN 2 log N + 
N 2 log kN log iV) expected time by recomputing the funnels each time the decision problem is computed. 
Note that 0(N 2 ) storage is required for the red-blue intersections algorithm (cf. Theorem [53]) . 
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7. Geodesic Hausdorff Distance 

Hausdorff distance is a similarity metric commonly used to compare sets of points or 
sets of line segments. The directed geodesic Hausdorff distance can be formally defined 
as 8h(A, B) = sup aG ^ inf&gs d(a, 6), where A and B are sets and d(a,b) is the geodesic 
distance between a and b (see [H [5]). The undirected geodesic Hausdorff distance is the 
larger of the two directed distances: 5jj(A, B) = m&x.(6H(A, B), 8h(B,A)). 

Theorem 7.1. 5h(A,B) for point sets A,B inside a simple polygon P can be computed in 
0((k + N) log (A; + N)) time and 0(k + N) space, where N is the larger of the complexities 
of A and B and k is the complexity of P. If A and B are sets of line segments, 5h(A,B) 
can be computed in 0{kN 2 a(kN)\ogkN) time and 0(kNa(kN) log kN) space. 

Proof Sketch. A geodesic Voronoi diagram [22] finds nearest neighbors when A and B are 
point sets. When A and B are sets of line segments, all nearest neighbors for a line segment 
can be found by computing a lower envelope [T| of O(N) hourglass distance functions. The 
largest nearest neighbor distance over all line segments is 5h(A,B). m 



8. Conclusion 

To compute the geodesic Frechet distance between two polygonal curves inside a simple 
polygon, we have proven that the free space inside a geodesic cell is x-monotone, y-monotone, 
and connected. By extending the shortest path algorithms of [131 HE], the boundaries of a 
single free space cell can be computed in logarithmic time, and this leads to an efficient 
algorithm for the geodesic Frechet decision problem. 

A randomized algorithm based on red-blue intersections solves the geodesic Frechet 
optimization problem in lieu of the standard parametric search approach. The randomized 
algorithm is also a practical alternative to parametric search for the non-geodesic Frechet 
distance in arbitrary dimensions. 

We can compute the geodesic Frechet distance between two polygonal curves A and B 
inside a simple bounding polygon P in 0(k+N 2 log kN log N) expected time, where N is the 
larger of the complexities of A and B and k is the complexity of P. In the expected case, the 
randomized optimization algorithm is an order of magnitude faster than a straightforward 
parametric search that uses Cole's [9] optimization to sort 0(kN 2 ) values. 

The geodesic Frechet distance in a polygonal domain with obstacles enforces a homotopy 
on the leash. It can be computed in the same manner as the geodesic Frechet distance inside 
a simple polygon after computing cell boundary funnels using Hershberger and Snoeyink's 
homotopic shortest paths algorithm [16J. Future work could attempt to compute these 
funnels in Oilogk) time instead of 0(k) time. The geodesic Hausdorff distance for point 
sets inside a simple polygon can be computed using geodesic Voronoi diagrams. The geodesic 
Hausdorff distance for line segments can be computed using lower envelopes; future work 
could speed up this algorithm by developing a geodesic Voronoi diagram for line segments. 
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